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From  examination  of  EKXL  code  run  output  and  of  analytic  studies  of  particle 
density  distributions,  a  picture  emerges  of  the  behavior  of  ions  and  electrons  in  flow 

fjj. 

through  simple  Polywell  systems. 

This  shows  that  the  region  within  the  momentum-limited  core  is  net  positively 
charged,  and  always  develops  a  virtual  anode.  If  the  ion  cij^ent,  for  given  electron  drive 
current,  is  allowed  to  become  sufficiently  large,  this  virtual  anode  increases  in  height  above 
the  potential  well  minimum  and  the  core  radius  broadens  accordingly.  This  broadening 
follows  the  simple  algorithm 
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where  E.^^  is  the  stable  transverse  ion  energy  at  the  outer  boundary  after  core,  edge,  and 
mantle  collision  processes,  and  E.^  is  the  potential  well  depth  at  the  center  of  the  core  (i.e. 
electron  maximum  injection  energy  minus  virtual  anode  height).  This  broadening  leads  to 
reduced  core  density,  varying  as  (1/r^),  for  given  drive  conditions. 

However,  since  the  ions  approaching  the  core  now  must  "climb"  the  increasing 
anode  "hill",  they  will  be  slowed  down,  and  their  radial  motion  will  decrease  with  energy 
given  up  to  the  central  virtual  anode  potential.  Thus,  their  convergence  will  be  more  rapid 
than  (the  usual)  geometric  inverse-square  (1/r  )  dependence  that  characterizes  flow 
throughout  most  of  the  device,  and  the  ion  density  will  increase  more  rapidly  than 
otherwise  expected.  This  increased  densification  tends  to  oflset  the  geometric  reduction  in 
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By  "limple"  syatenu  ii  meant  potential  welle  whose  shape  is  describable  by  a  simple  power 
law  expression  of  the  form  <r>™;  with  m  as  m  s  3  for  the  Polywell^'"/SCIF  systems  of  interest  here 
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density  due  to  increased  core  radius.  The  net  result  is  that  fusion  reaction  rates  vary  only 
slowly  with  increasing  anode  height  up  to  about  half  the  wdl  depth.  Beyond  this  point  the 
reaction  rates  steadily  decrease  due  to  decreasing  fusion  reaction  cross-sections  with  the 
decreasing  ion  kinetic  energy.  This  feature  of  ion  flow  is  discussed  in  more  detail  in  a 
forthcoming  EMC2  Technical  Note.^ 


In  the  region  outside  the  core  the  plasma  is  near-neutral  (though  still  slightly 
positive  out  to  the  well  minimum,  at  r  =:  r^)  and  the  ion  and  electron  density  follow  each 
other  very  closely  (i.e.  within  <  lE-5  of  local  density,  for  all  densities  of  interest), 
decreasing  about  as  the  inverse-square  of  the  increasing  radius.  The  variation  from 
neutrality  is  such  that  the  net  charge  density  varies  approximately  linearly  with  radius, 
decreasing  to  zero  at  r^.  Beyond  r^  the  system  is  nearly  neutral,  but  increasingly  slightly 
negative,  out  to  a  radius  r^  at  which  the  total  energy  of  the  average  electron  is  found  to  be 
in  well  potential  energy  and  transverse  kinetic  motion  (i.e.  where  the  average  electron  has 
no  radial  kinetic  energy).  This  has  been  called  the  "stagnation  radius"  for  electron  flow;^ 
in  principle,  in  the  absence  of  ions,  the  average  electron  would  never  reach  a  radius  smaller 
than  r^  This  radius  is  just  the  mean  convergence  radius  of  electrons  with  transverse  energy 
spread  ±  injection  energy  E^,  and  is  given  by 
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For  example,  if  the  electrons  are  given  transverse  energy  at  the  outer  boundary  of  Por  /  I 

E _ =  0.6  E^,  then  the  value  of  <z^  =  (r^/R)  will  be  <rj>  =  0.775.  For  a  diffuse:  ^ 
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reflection  distribution  of  electrons  from  the  outer  boundary  region,  E  =  E  /3  and! 

°  ’  ep«rp  o'  j 

<rj>  =  0.578. 
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The  gross  density  of  both  species  continues  to  drop  approximately  as  the  inverse 
square  from  r^  out  to  a  radius  r  =  at  which  both  the  ion  and  electron  density  begin  to 
increase  rapidly  as  the  ions  lose  energy  to  the  potential  well  and  the  electrons  pile  up  to 
form  the  main  body  of  this  potential  well.  This  region  from  r^  <  r  <  R  is  dominated  by 
electron  motion.  Outside  rj^  the  electron  density  increases  much  faster  with  radius  than 
does  the  ion  density,  and  the  system  becomes  highly  net  neg^ive  thus  ensuring  the 
establishment  of  the  deep  negative  well  required  for  ion  acceleration. 


Ion  motion  dominates  the  flow  within  r  <  Between  r^  <  r  <  rj^  the  transition 
changes  from  ion-driven  flow  (at  r  <  r^)  to  electron— driven  flow  (at  r  >  r^).  The  density 

A 

at  rj^  is  found  from  the  distributions  derived  by  Krall  to  be 


(3) 


for  a  simple  power-law  well  with  m  =  3,  where  the  potential  distribution  is  given  by 
e0(r)  =  e^^[l-(r/R)^]  =  a^E^[l-(r/R)^],  defining  the  well  depth  factor 

This  work  gave  ion  density  distributions  in  three  regions  across  the  system.  These 
were  in  the  "core”  r  <  r^,  the  "mantle"  r^  <  r  <  rj^,  and  the  outer  region  r  >  r^^.  It  also 
showed  that  the  ion  density  within  a  very  small  radial  distance  in  from  the  "edge"  at  r  =  R 
would  have  to  rise  slightly  in  order  to  preserve  flow  and  current  continuity  within  the 
electron-driven  potential  well.  The  "true"  monotonic  potential  well  thus  starts  at  a  radial 
position  r^  slightly  less  than  R;  <<  r^  <a  R.  From  this  body  of  work  it  is  possible  to 
obtain  the  properties  of  the  flow  at  the  several  critical  radii  just  discussed. 
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Consider  the  "edge"  region,  and  suppose  that  the  potential  well  shape  is  a  simple 
(m)  power-law  here,  so  that  ion  and  electron  energies  would  be  given  by 


Ej(r)  =  (l/2)m.y  *  =  E,(l  -  <r>“) 

(4a) 

E.(0  =  (1/2K =  E^(<r>") 

(4b) 

where  <r>  =  (r/R).  The  radius  <rm>  at  which  the  ion  density  will  be  a  maximum  is  that 
at  which  the  potential  well  has  fallen  &om  r  =:  R  to  a  depth  equal  to  the  initial  spread  in 
energy  of  the  ion  input,  here  assumed  to  be  uniformly  distributed  &om  transverse  to 
radial.  The  convergence  radius  is  determined  by  this  spread  as  <r^>  =  (E.^/E^)®  ®,  so  that 
the  "edge"  radius  is  given  from  the  ion  eq.  (4a)  as 
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which  is  very  small  for  all  convergence  ratios  and  power  law  exponents  of  interest.  For 
example,  for  m  =  3  and  <r^>  =  lEl-2,  eq.  (6)  gives  =  1  -  3.3E-5. 


Now  consider  the  radius  <r^>  which  electron  and  ion  velocities  would  be  equal, 
assuming  the  validity  of  the  radial  flow  model  used  here.  This  is  readily  found  from  eqs. 
(4)  to  be  simply 
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For  ions  of  m  =  3  and  mass  2m^  (e.g.  D)  this  gives  <r^>  =  0.0650.  Of  course,  in  any  real 
system  the  existence  of  finite  transverse  momenta  and  their  conservation  will  render  this 
particular  numerical  result  invalid,  but  it  is  indicative  of  the  nature  of  the  ion  and  electron 
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flow  in  the  system. 


In  considerable  contrast  is  the  radius  <r  >  at  which  ion  and  electron  momenta  are 

p 

equal.  Again,  from  eqs.  (4)  this  is  found  to  be 


<rp>  =  l~(l/m) 


(7) 


which  yields  <rp>  =  1  -  9.07Bi-S  for  D  ions.  Note  that  this  is  very  near  to  the  value  for 
the  radius  of  maximum  ion  density,  above.  Finally,  solving  eqs.  (4)  for  the  radius 
<r^>  at  which  ion  and  electron  energies  are  equal  gives 

s  * 

<r^>  =  0.794  (8) 


The  critical  radius  <r|^>  at  which  the  ion  density  reaches  a  minimum  is  found  by 
differentiating  the  equation  (given  previously  )  for  ion  density  in  the  region  between  <rj^> 
and  <r^>-  This  yields  the  result  that 

<r^>  =  0.830  (10) 

not  far  different  from  the  value  for  <r^>  just  obtained. 

The  ion  densities  expected  from  this  simple  model  at 
summarized  below  in  terms  of  Krall’s  edge  density  n^^,  and 


these  critical  radii  are 
the  core  density  n^. 
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ION  DENSITY 

IN  TERMS  OF  nj^ 

IN  TERMS  OF  n 

c 

At  core  edge 

n.  =  (1.5/<r^>)nj^ 

n.  =  n 

1  c 

for  actual  density  distribution  in  system 

At  radius  of  min¬ 

n.  =  (3.327<r^>)nj^ 

EL  =  (2.218<r^>*)n^ 

imum  density  <rj^>, 

<r,  >  =  0.830 

k 

2  ^ 

for  (1/r  )  density  distribution  at  r  <  r^ 

n.  =  (2.180<r^>)njj^ 

n.  =  (1.453<r^>^)n^ 

At  radius  of  max¬ 
imum  density  <r  >, 
m 

=  (l-<re>V3) 

n.  =  (1.31)nj^ 

n.  =  (0.873<r^>)n^ 

At  system  boundary 

II 

n.  =  (0.667<r^>)n^ 

These  formulae  ate  of  some  interest,  as  they  illustrate  the  general  nature  of 
variation  of  ion  density  across  the  device  when  it  is  operating  with  a  stable,  filled  potential 
well  of  simple  shape.  In  particular  it  is  interesting  to  note  that  the  core  ion  density  is 
higher  than  the  density  at  injection  by  a  factor  of  the  order  of  (R/r^),  and  higher  than  the 

9 

minimum  density  in  the  system  by  (R/r^)  .  The  electron  density  follows  the  ion  density 
closely  for  all  radii  r  <  r^.  Beyond  r^^  the  electron  density  increases  markedly,  becoming 
comparable  to  the  core  density  at  the  system  boundary. 

In  actual  fact,  this  simple  model  becomes  progressively  more  in  error  as  the  outer 
boundary  is  approached,  because  the  B  field  -  which  determines  the  potential  shape  —  does 
not  follow  a  simple  power  law  formula.  Rather  it  varies  about  as 
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B  =  B. 


2<r> 


m 


(n) 


This  "roUoTer"  field  shape  leads  to  aflat  B  field  gradient  at  the  system  edge,  which 
strongly  affects  the  density  distribution  in  this  repon.  The  effect  is  to  flatten  the  electron 
distribution  considerably  from  its  sharply  rising  character  under  the  simple  model  used 
here.  This  means  that  the  ion  density  will  also  not  rise  as  rapidly  with  increasing  outer 
radius  position,  and  that  aU  calculations  whose  answers  are  driren  by  the  high  edge 
densities  of  the  simple  model  are  suspect  and  probably  in  error.  Within  r  <  r^  both  models 
are  fairly  close  and  the  distributions  given  by  the  simple  model  here  may  be  tahen  as 
reasonable  approximations  to  real  effects.  It  is  only  in  the  outermost  regions  of  the  system 
that  significant  complexities  arise;  and  these  will  not  be  resolved  short  of  good  experiments 
that  show  particle  and  field  distributions. 
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